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The annihilation of positronium is discussed as a means of obtaining a precise measurement of 
the Compton wavelength of the electron. In particular the spectral line shape of the annihilation radia- 
tion is calculated for positronium thermalized in helium gas cooled to the region of its critical point 
with an applied magnetic field sufficient to yield a three photon to two photon branching ratio of unity 
for the 3 S component. Under such conditions the rapid thermalization of the positronium and the long 
half-life of two photon decay out of the predominantly triplet state make possible the production of 
rather narrow spectral distributions. An exact solution of the Zeeman effect in positronium is included. 
This study suggests that a wavelength determination beyond the one part per million level is attainable 
within the range of present technology. 
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1. Introduction 

X-ray data have been excluded from recent adjust- 
ments of the atomic constants [1J 1 since at their then 
current precision they would have carried typically 
about one fifth the weight assigned to other measure- 
ments. Simply improving the precision of the x-ray 
data by -a factor of three would lead to equal weight. 
Recent improvements in the internal consistency of 
the x-ray scale by Bearden et al. [2], permit a further 
sharpening by x-ray determinations of the constants to 
limits set by the uncertainty in the conversion factor 
between the x-ray wavelength scale and the denned 
length. Were these potentialities realized today, then 
the x-ray data would be significant in an adjustment of 
the fine structure constant and the proton moment [3J. 

The precision of x-ray measurements of the atomic 
constants is -limited by the precision with which it is 
possible to define an x-ray wavelength scale. The scale, 
of course, depends on the precision with which it is 
possible to establish the wavelength of an individual 
x-ray line. Ordinary x-ray lines have full widths at 
half maximum upwards of some 300 ppm (parts per 
million). We have no theoretical model for their shape 
and hence little hope of specifying the wavelength 
of any of their features (peak, centroid, or median) 
to better than 1 ppm. This problem has been studied 
by J. S. Thomsen and F. Y. Yap [4). 
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Beside the traditional x-ray determinations of the 
constants there are phenomena which either occur 
spontaneously or can be provoked to occur which 
are considerably sharper than natural x-ray lines. 
For example, there are nuclear transitions which have 
quite narrow profiles [5] and there are crystal lattices 
of extraordinary perfection [6] which occur as, so to 
speak, spontaneous phenomena. Other phenomena 
which are also quite sharp can be produced by ex- 
perimental design. One of these, namely radiation 
from the annihilation of positronium as narrowed by 
thermalization in a magnetic field is the subject of 
the present paper. 

2. Positron-Electron Annihilation 

If a positron and an electron initially at rest an- 
nihilate with the emission of two photons, then by 
conservation of energy, 



2^-=2m c\ 



(1) 



since conservation of linear momentum implies that 
the wavelengths of the two photons are equal. Solving 
eq (1) for A. yields the Compton wavelength: 



k A = 



h 

moc 



(2) 



where m () is the rest mass of either the electron or 
the positron. 
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If, however, the positron-electron system has non- 
zero momentum at the instant of annihilation, the two 
photons in general will not have the same wavelength, 
so that spectroscopic observation of the radiation 
shows a broadened annihilation line. Examination of 
this process shows that the line widths resulting from 
the annihilation of free positrons with electrons in 
matter are much too great (typically in excess of 1000 
ppm) to allow a precision (i.e., approx. 1 ppm) measure- 
ment. 

If a flux of positrons is incident upon certain ma- 
terials, a significant fraction of these incoming par- 
ticles will combine with electrons in the substance to 
form hydrogen-like atoms of positronium [7]. The peak 
of the annihilation line is shifted because the right- 
hand side of eq (1) is reduced by the binding energy 
of the positronium atom, while the width of the line 
is determined by the momentum of the positronium 
atom at the time of annihilation. 

Positronium exists in either singlet (para-) or triplet 
(ortho-) states; the lifetime against annihilation of 
ortho-positronium is two orders of magnitude greater 
than that of para-positronium. The triplet lifetime is 
sufficiently great to allow an ortho-positronium atom 
to reach thermal equilibrium with its surroundings 
before annihilating. If this thermalization process is 
carried out in the presence of an external magnetic 
field, then an appreciable fraction of the long-lived 
positronium component will decay via the two photon 
channel because the field causes mixing of the singlet 
and triplet states. If the temperature of the surround- 
ings is sufficiently low, then the resulting annihi- 
lation line will be narrow. 

It is the purpose of this paper to examine in detail 
the line shape of the radiation produced in this process. 
In particular we shall consider the lifetime of the two 
photon annihilation of positronium in the presence of a 
magnetic field; combining this with appropriate ex- 
pressions for the thermalization of the positronium 
and for the thermal distribution function of the sur- 
roundings yields an expression for the line shape. 



can represent the stationary ground state as an eigen- 
state 



S&oUrm — EfUftt 



(3) 



where F = 0, 1 is the quantum number denoting the 
spin state of the atom, and m is the quantum number 
associated with z-component of the spin angular mo- 
mentum. The states Ufm are eigenfunctions of the spin 
operators (S e + S p ) 2 and {S ez + S pz ), where the sub- 
scripts e and p denote the electron and positron, re- 
spectively. These states, written in terms of the spin 
states of the individual particles, are 



aii = -4WT')Mt) 



(4a) 



^o 



hio = —= M t )M 1 ) + M i K( t )] (4b) 

V2 



Ui-i =^ Ue( i)Up( 1) 
^o 



(4c) 



and "oo = —= [u e ( t )iip( I ) - M I )u p ( f )] (4d) 

V2 



where ^o is the hydrogenic ground state wave function 
which results if spin is neglected. The hfs separation 
has been evaluated to be 

A£ = £i -£ = 8xlO- 4 eV. 

If an external magnetic field H is applied to the 
positronium atom, then the Hamiltonian becomes 



& = &0 + &H 



(5) 



where, in terms of magnetic moment operators /x for 
the positron and electron, 



$Qh=— (#*e+Mp)'H. 



(6) 



3. Zeeman Effect of Positronium 

The hyperfine separation between the singlet and 
triplet energy levels of positronium in its ground state 
is much greater than the corresponding splitting in the 
hydrogen atom for two reasons [8]. In the first place 
the dipole-dipole interaction between the electron and 
positron is proportional to the magnetic moment of 
the positron, and therefore is increased over the inter- 
action in hydrogen by the ratio of the Bohr magneton 
to the nuclear magneton. Furthermore, the hyperfine 
splitting is increased by an amount that is approxi- 
mately equal in magnitude to the dipole-dipole term 
by quantum electrodynamic effects involving virtual 
annihilation of the pair. 

If §o represents the Hamiltonian (including the 
hyperfine interaction) of the positronium, then one 



Since the magnetic moments of the electron and 
positron are equal (approximately) in magnitude to 
the Bohr magneton but opposite in direction with 
respect to the spin vectors, and since their g factors 
are equal, then S*)h can be written in terms of Pauli 
spin matrices cr z as 



_ fiogH 

&H — (o-ez — <Tpz). 



(7) 



Some treatments [9] of the Zeeman effect in posi- 
tronium have employed perturbation theory; however, 
an exact solution is readily obtained by well-known 
methods [10]. 

One can write down the Hamiltonian matrix by 
employing the fact that a z u( / [)=u(\) and a z u( I ) 
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= — u( I); the operation of &h in eq (7) on the wave 
functions in eq (4a-d) yields 



$Qh\uu =!Qh\ui-i > = 

$h\u\0 > =IJLogH\uoo > 

&h\uoq > = [M>gH\u l0 >. (8) 



and 



As a result, the Hamiltonian matrix including the mag- 
netic interaction is 



<$> = 







fiogH 




(9) 



The secular equation for (,<p) is 

(£,-A) 2 {(£.-A)(£o-A)-(|Li (# //) 2 }=0 

with roots 

1 



where 



and 



k^EuEuE^-AEVl + x* 



E = |(£, + £ ) 



(10) 



(11) 



IM)gH 

AEI2' 



Thus the diagonalized matrix is 
'Ei 

E + ^Vl+x 2 



(^)diag-l 
















o E-^Vi-i-x 2 



(12) 



Note that the degeneracy between the triplet m = ±l 
states is not removed by #//. Furthermore, the energy 
of these states is not affected by the field, since these 
states have zero net magnetic moment because the 
positron and electron are oppositely charged particles. 
The lack of a linear Zeeman effect, and the existence 
of a quadratic effect in the m = state can be described 
qualitatively as follows. The triplet m = state as 
mentioned above has no magnetic moment and is 
not affected in first order. The singlet state, though 
it possesses a magnetic moment, is a zero spin state; 



since there is no axis of quantization, the atom 
"tumbles" in the field. Both states do exhibit an 
energy shift because the applied field induces a 
preferential alinement of the individual particle 
magnetic moments as is indicated by the H 2 variation 
(in the limit, this induced alinement becomes the 
Paschen-Back effect). 

The wave functions for the perturbed positronium 
atom can now be determined by carrying out the matrix 
multiplication 



<S><£> = <£>diaK<S> 



(13) 



where (S) is the transformation matrix that diag- 
onalizes (£>), and employing the fact that (S) must be 
unitary. This operation leads to the following wave 
functions: 




^here 



«-4:(l+- l 



\fi V vTT 



J 



and 



V2 



VT 



1/2 



1/2 



(15a) 



(15b) 



Equation (14) shows that the application of a mag- 
netic field mixes the states Uio and uoo- Thus \jho and 
i/ibo each decay by both the two- and three-photon 
channels. 

If the annihilation rate of the state Uio is yr and that 
for the state uoo is y s , then the mean annihilation rate 
for the mixed state 1//10 is 



yio = a 2 y T + b 2 y s . 
Substitution of eq (15) into the above yields 

yr+ys\ , l / yr — y> 



yio = 



+ ■ 



Vl^K 



(16) 



(17) 



The branching ratio for the 1//10 state (i.e., the ratio 
of the number of two photon decays to three photon 
decays) is 



B 



_(b\*ys _ (2 + x 2 )-2\/l+x 2 y s 



-i 



y T 



Jt 



(18) 
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In figure 1 are plotted eqs (17) and (18) as functions 
of the applied field H up to 2177 G, at which point 
Z?io equals unity. 
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FIGURE 1. Quenching of triplet positronium at low fields. 



FIGURE 2. (a) Two-photon annihilation of positronium in motion, 
(b) Momentum triangle for the process (a). 



Line Profile for Positronium in 
energetic Moderator 



a Mono- 



In this section we shall consider the case in which a 
flux of positronium atoms is incident in a moderating 
medium possessing a single thermal energy. We must 
first relate the energy separation of the two photons 
to their emission directions. Figure 2a is a represen- 
tation of the annihilation of a positron atom with linear 
momentum p into photons hv\ and hv 2 which are 
emitted at an angle e with respect to each other. The 
photon hv\ is emitted at an angle with respect to the 
direction of the incident positronium atom. From the 
momentum triangle of figure 2b one can write 



pI-p\- 



p 2 — 2pp\ cos 6. 



(19) 



and subsequently substituting eqs (20) into the result 
leads to 



8 = - 



pc 

2 



pc_ (l + 28 



E„ 



£o 



COS 



(21) 



If the positronium atom is nearly in thermal equilibrium 
with its surroundings at annihilation, then 



pc 
E 



< 1 



and 






Letting hv\ be the photon of greater energy, then one 
can write 



as a result, eq (21) can be written 



and 



hvi=lE + d 



hv> = | £0—8, 



(20a) 



(20b) 



where £0 is the total energy of the positronium atom 
prior to annihilation. Thus multiplying eq (19) by c 2 



pc 



cos = 8 max cos 9 . 



(22) 



It has been shown by Heitler [11] that the emission 
of annihilation photons is spherically symmetric in 
the limit p=0; thus the distribution will be very nearly 
isotropic in the case at hand — that of thermalized 
positronium. The intensity distribution as a function 
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of 6 will thus be proportional to the solid angle sub- 
tended at the annihilation site: 3 



1 



IedO=^sinOdO, Q<\ 



(23) 



Then the intensity distribution for 8 can be written 



dft 
hd8 = I 9 dO=Io—d8,8<8 n 

do 



(24) 



The quantity d6/d8 can be evaluated from eq (22) so 
that the distribution function I§d8 in eq (24) becomes 



hd8 — pc, 8 < 6 max . 



(25) 



Thus the annihilation of a positronium atom traveling 
with momentum p is equally likely to produce photons 
in the range from £"0/2 up to Eq[2 + pc/2. 

It is now possible to determine the annihilation line 
shape from a positronium flux in a monoenergetic 
moderator. The number of Ps atoms annihilating during 
the time interval (t, t + dt) is 



Ntdt = yioN e-v»> f dt. 



(26) 



Annihilations during this time interval contribute an 
increase of intensity dh which is uniform for all 5 
from up to 8 max , according to eq (25). Since the area 
under the distribution curve is equal to the number of 
atoms annihilating, then 



N t dt = 8 max dl 8 . 



(27) 



We can determine the total intensity distribution by 
adding up all the diss from £ = up to t$ such that 

O Ornax* 



h 



[ h Ntdt 
Jo P° 



(28) 



The momentum p= v2mE is a function of the time; 
this thermalization function has been found [12] to be 



E = E m coth 2 (/3 + a*) 



(29) 



where E m is the moderator energy, f3 depends upon the 
energy of the positronium atoms at t = 0, and a de- 
pends upon the properties of the moderator. Defining 



-fi + at 



eq (28) can be written as 

Jo 



d tanh udu. 



(30) 



1 In reality 6 does not quite reach tt/2; this is merely an indication of the approximation 
of spherical symmetry, and will not affect the results significantly. 



This integral can be approximated in the case for which 
yio/o: is small. The details of this calculation appear 
in the appendix. The resulting distribution function is 



:/o[l 



1 



= /o 



8>8 T 



0<8<8 T 



(31) 



where 8r is the maximum energy shift for the annihi- 
lation of a completely thermalized positronium atom. 
In the case of two-photon annihilation of the i/>io 
component of positronium using helium at its critical 
point (7V = 5.2°K, /; = 2.3atm) as a moderator and 
an applied magnetic field of 2177 G (so that Z?i =l), 
the quantity yio/a is quite small: 



yio 

a 



Hecrit = 8.5X10-3. 



As a result the function /g (eq (31)) falls very rapidly 
for 8 > 6Y, so that this leads to a distribution function 
that is nearly rectangular; viz, at 8/8r=l,/=/o, 
while at 8/8t= 1.001, I=h X 10~ 3 . 

The physical interpretation of this result is fairly 
simple. Since y U) is the two photon annihilation rate 
of the long-lived component of positronium, then 
yfo is the mean lifetime of the atom; similarly, a~ ] 
represents the characteristic time required for the 
Ps atom to thermalize with its surroundings. Thus the 
ratio yio/a is just the ratio of thermalization time to 
mean life. When this ratio is small, then thermalization 
usually is completed long before the positronium 
atom decays. Since the calculations of this section 
concern the somewhat artificial case of a mono- 
energetic moderator, the resulting rectangular, dis- 
tribution in the preceding paragraph indicates that 
the vast majority of the positronium atoms decay after 
they have become fully thermalized, so that the re- 
sulting intensity distribution is simply the charac- 
teristic rectangular distribution arising from mono- 
energetic positronium. 



5. Line Profile for Thermalized Positronium 

Let us now consider annihilation in a moderator 
at temperature T characterized bv a Maxwell dis- 
tribution of velocities: 



N v dv = Av 2 e ^kT dv. 



(32) 



If the thermalization time a' 1 is much less than the 
mean life yfo 1 . then most Ps atoms are in thermal 
equilibrium with the moderator at the time of anni- 
hilation. As a result there will be a Maxwell distribu- 
tion of positronium atoms; atoms with velocities 
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in the interval (v, v+dv) will produce a rectangular life time sufficiently great to fulfill the conditions 

for the Gaussian line shape described above. It would 
thus appear possible to carry out an annihilation 



distribution according to eq (31), so that 
ddI=N v dv 



133) 



where 8 is the width of the distribution and is the 
maximum energy shift for photons produced in the 
decay of a Ps atom moving with velocity v. If one 
makes the substitution 8—1/2 mvc, eq (33) can be 
integrated to yield the intensity distribution of anni- 
hilation photons: 



experiment in which the line width is sufficiently 
narrow to allow a 1 ppm determination of \a. 

7. Appendix 

In order to evaluate the integral of eq (30), 



-I 






2kT dv 



-/■ 



tanh udu, 



(A-l) 



8/mc m,VC 

= Ke 



2kT\mc) 



(34) 



for the case in which € = (yiol<x) < 1, we make the 
substitution 



The quantity 8 can be expressed in terms of the posi- 
tronium kinetic energy as 



77 — In cosh u. 



(A-2) 



8 = cV2mE. 



(35) 



Denoting as 8r that value of 8 corresponding to the 
mean thermal energy at T, then one can write that 



dr — 



yJ2m(^kT^ = yl(^kTJE . (36) 



Since £o °" 2moc 2 , then eq (34) becomes 
I = Ke 



UfcJ 



(37) 



The line profile resulting from the annihilation process 
is a Gaussian with a half width at half intensity of 



The expression tanh udu can be expressed in terms 
of 7): 

, , d(cosh u) in 1X7. 

tanh udu = r =d(ln cosh u) = dr\. (A-3) 

From (A-2) we can write 

e 17 = cosh u 
which leads to 

u — cosh -1 ^ 
and therefore (A-l) becomes 

S— I e _ecosh " le " dr\. 



(A-4) 



(A-5) 



(A-6) 



81/2 — 8t 



Employing an identity for cosh" l yields 



In 2 = 0.68 8r. 



(38) 



6. Summary 



cosh- 1 e r > = In (e 17 -fVe 2 ^- 1 ) 

= 7} + In (l + Vl-e- 2r 



and therefore 

e -e cosh-i^ = ^-^d + Vl-e" 2 ^)- 



Since e is small, eq (A-8) can be expanded: 
(l 4- Vl-e- 2r >)- e =l-€ In (l -Vl-e- 2r 

+ |"ln 2 (l + Vl-e- 2 ^)-. . . . 



It has been demonstrated that if a moderator is 
such that the thermalization time of positronium is 
much less than the mean life against annihilation, 
the two photon annihilation of positronium in this 
moderator at temperature T results in an annihila- 
tion line that is a Gaussian, the width of which depends 
upon T. These conditions are fulfilled for helium 
near its critical point for positronium lifetimes of the 
order of magnitude of the triplet state lifetime. Since 
triplet positronium decays with the emission of 
three photons, whereas a determination of the Compton But note that 
wavelength is feasible only in the case of two photon 
decay, one must carry out such an experiment in } n (l + Vl — e~ 2 v 
the presence of a magnetic field so that the mixing 
of singlet and triplet states give rise to a long-lived 
two photon process. It has been shown that for fields and that one can write 
of approximately 2 kG, one half of the long-lived 
decays will occur via the two photon channel with a (l -f Vl — e~ 2v )~ e = 1 
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(A-7) 



(A-8) 






-2t)\*/2 



<t>€ 



(A-9) 



(A-10) 



(A-ll) 



in w 


hich 








e^ eln (l+Vl — e" 


u) 


Now 


the ] 


ntegral (A— 6) becomes 




9 


= e-^l-f *€)A7 = 


-( 


where 







,e(l — e- 2r >yi 2 . 



-4>€) 



(A-12) 



•-/ 



<be 



Rewriting (A- 
tion 3. 



[1] E. R. Cohen and J. W. 
537 (1965). 



drj / e-^dr), 0< <&< 1. 
12) in terms of u yields eq (31) in sec- 
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